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Abstract 

We study several issues related to the use of effective field theories in QCD at large baryon 
density. We show that the power counting is complicated by the appearance of two scales inside 
loop integrals. Hard dense loops involve the large scale and lead to phenomena such as screening 
and damping at the scale g^. Soft loops only involve small scales and lead to superfluidity and 
non-Fermi liquid behavior at exponentially small scales. Four-fermion operators in the effective 
theory are suppressed by powers of but they get enhanced by hard loops. As a consequence 
their contribution to the pairing gap is only suppressed by powers of the coupling constant, and 
not powers of l//i. We determine the coefficients of four-fermion operators in the effective theory 
by matching quark-quark scattering amplitudes. Finally, we introduce a perturbative scheme for 
computing corrections to the gap parameter in the superfluid phase. 



I. INTRODUCTION 



The study of baryonic matter in the regime of high baryon density has led to the discovery 
of new phases of strongly interacting_matterj_ su'^i. color superconducting quark matter 
and color-flavor locked matter [l|, O, ISl |6|, |7, 8|. These phases are not only relevant to 




the structure of compact astrophysical objects [9|, llO|, but they also provide a theoretical 
laboratory in which complicated QCD phenomena, such as chiral symmetry breaking and 
the formation of a mass gap, can be studied in a weakly coupled setting jn|. In order 
to exploit these opportunities we would like to develop a systematic framework that will 
allow us to determine the exact nature of the phase diagram as a function of the density, 
temperature, the quark masses, and the lepton chemical potentials, and to compute the low 
energy properties of these phases. 

If the density is large then the Fermi momentum is much bigger than the QCD scale, 
p_F = ^ ^QCDi and it would seem that such a framework is provided by perturbative 
QCD. It is clear, however, that a naive expansion in powers of as is not sufficient. First 
of all, it is well known that resummation is required in order to make the perturbative 



expansion in a many body system well defined 12] • It is also well known that additional 
problems arise in systems with unscreened transverse gauge boson interactions 
a degenerate Fermi system the effect of the BCS or other pairing instabilities have to be 
taken into account. And finally, in systems with broken global symmetries, the low energy 
properties of the system are governed by collective modes that carry the quantum numbers 
of the broken generators. 

In order to address these problems it is natural to exploit the separation of scales provided 
by /i ^ ^ Aqcd in the normal phase, or fi ^ g^i ^ A ^ Aqcd in the superfiuid phase. 
An effective field theory approach to phenomena near the Fermi surface was suggested by 



long 



This approach was applied to the gap equation [l^l, Goldstone boson masses 



gluon dispersion relations |20|, and a number of other problems. For a review and 
further references, see j^]. Even though a number of interesting results have been obtained 
there are a number of important conceptual issues that are not very well understood. These 
issues concern power counting, renormalization and matching. In this work we would like to 
study some of these issues in more detail. This paper is organized as follows. In sectionlTTlwe 
give a review of the standard hard thermal loop approximation applied to dense matter. In 
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Section ITTTl we introduce the high density and effective theory (HDET) and in Section HVl we 
show how hard dense loops (HDLs) arise in the effective theory. In Section |3 we show that 
the effective theory contains a new class of diagrams which we shall call soft dense loops. In 
Sections I VII and IVIII we study the gap equation at leading and next-to- leading order and in 
Section IV nil we study matching conditions for four-fermion operators. 



II. HARD DENSE LOOPS 



In this section we review the calculation of hard dense loop contributions to QCD Greens 



functions 
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2J|. The hard dense loop limit corresponds to soft external momenta, 
\p\ <^ fi. In this limit the main medium contribution comes from hard loop momenta, 
\k\ ~ fi. Our main purpose in this section is to present a simple rederivation of the main 
results that will allow us to compare the hard dense loop approximation with the high 
density effective theory which is discussed in section IIIII 

Hard dense loops can be calculated most simply by writing the free fermion propagator 
at 7^ in the form 





where = ±(|A;| =l= /i) and K = (1, k), K = (1, —k) with k = k/\k\. The two poles of the 
propagator in equ. (jT]) correspond to particle and anti-particle states. Note that = -^7o/2 
and A^ = are projection operators on positive and negative energy solutions of the 

free Dirac equation. Also observe that (-^)^ = (-^)^ = which is very useful in evaluating 
Dirac traces in the hard dense loop approximation. In our convention the energy is measured 
relative to the Fermi energy fi. In order to derive HDL Greens functions it is useful to define 
the energy of a fermion as u = ko + fi. The HDL limit corresponds to uj,\k\ <C n, whereas 
quasi-particles in the vicinity of the Fermi surface satisfy cj, ~ yU. 

As an example, let us consider the calculation of the gluon polarization function. The 
one-loop contribution is given by 

KiP) = Nf9Y J -^J^Mkh^S{k-p)). (2) 

If the propagator is written in terms of particle and anti-particle contributions as in equ. (|3ip 
we find three contributions to the polarization function, the particle-hole, particle-anti- 
particle, and anti-particle-anti-hole terms. The particle-hole contribution shown in Fig^ 
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UfMl.2,nH-l§A-J<M-J-^\ (3) 



p-K 

where = Nfg'^fi'^/lAir'^). This contribution is dominated by loop momenta I 
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25|. The 



close to the Fermi surface. However, the particle-hole term is not transverse 
particle-anti-particle term is 

Uf^ip) = 2m'5^' J ^ {5,0^.0 - K,K} ■ (4) 

Note that this term receives contributions from particles with momenta in the range < 
\k\ < /i, and not only particles in the vicinity of the Fermi surface. Also note that in the 
hard dense loop limit this term is just a constant, with no momentum dependence. The sum 
of the two contributions equ. and (jH) is given by 

n- (p) = 2m^5- / ^ Lo5.o - . (5) 



47r i '^^ p.K 

The polarization tensor equ. El) is transverse and agrees with the well known hard dense 
loop result. 

The fermion self energy in the hard dense loop approximation is given by 

with mj = Cpg'^j^'^/{8iT'^) and Cp = {N^ — l)/{2Nc). This result arises solely from the 
particle term in the fermion propagator equ. (|31|). The loop integral receives contributions 
from all momenta < |A;| < fi. It is important to emphasize that the HDL approximation 
corresponds to external momenta that are soft, \p\ <^ fi, whereas low energy fermionic 
excitations have momenta \p\ ~ fi. The use of the HDL fermion self energy is not reliable 
for low energy modes. It is nevertheless instructive to study the behavior of the HDL self 
energy near the Fermi surface. Using equ. © we find 



with 



= - Sv(p)7 ■ P (7) 



We observe that both Us and Sy are logarithmically divergent near the Fermi surface, but 
the quasi-particle dispersion relation and the wave function renormalization factor are well 
behaved. Indeed, for the particle mode 

s{p) = {So\p) - mr' ^ + • • • , (10) 

PO-UJ+ 2 

with 

.2 ^2 



.,(p)^|p1 + ^. Z,(p)^l + ^^l-log^^jj. (11) 

At one-loop order the quark-gluon vertex receives contributions from the abelian diagram 
Fig. and the non-abelian diagram Fig. Eb- In the HDL limit these two diagrams are 
identical up to a color factor. We find 

where Pi,P2 are the momenta of the two fermions and we have factored out the color structure 
X°'/2. The vertex correction is similar to the fermion self energy in the sense that the loop 
integral receives contribution from momenta < \k\ < fi. Also, we can study the effect 
of the HDL effective vertex for quasi-particles in the vicinity of the Fermi surface. Since 
both the incoming and outgoing Fermions are hard the interesting regime corresponds to 
soft gluon momenta. In this limit the matrix element of the free vertex between projectors 
on quasi-particle states with momentum p is (7^) = P^. For the HDL vertex correction we 
find 

/i J Att p' . K 

which has a logarithmic divergence near the Fermi surface. 

The gluonic three-point function in the HDL limit only receives contributions from the 
fermion loop diagram shown in Fig. Ek- We find 

rt,J,,,,r)..,r^2^'l§K,K^—JI^^ (14) 

The three-point function, as well as higher n-point functions, is dominated by momenta near 
the Fermi surface. Also, as opposed to the case of the two-point function, anti-particles only 
make sub-leading contributions. Higher n-point functions can be computed directly using 
the methods discussed here. Alternatively, as is the case for T 7^ 0, n > 3 point functions 
can be reconstructed using Ward identities or non-local effective actions. To summarize this 



section, we showed that gluon n-point functions in the HDL hmit are dominated by modes 
near the Fermi surface. The only exception is the two-point function which requires a contact 
term determined by modes with momenta between zero and the Fermi momentum. Fermion 
HDLs, on the other hand, are always determined by momenta in the range < \k\ < fi. We 
emphasized, however, that the HDL limit is not relevant for low energy fermionic excitations. 



III. HIGH DENSITY EFFECTIVE THEORY (HDET) 

In this work we wish to compare the hard dense loop results with the Greens functions 
obtained from an effective theory which describes low energy excitations in the vicinity of 
the Fermi surface. The leading order terms in the effective theory are 
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161, 



|18| 



C = Y.ijl{zv D)^, - \gI,GI, + . . . , (15) 

where = (1,-u). In the vicinity of the Fermi surface the relevant degrees of freedom are 
particle and hole excitations which move with the Fermi velocity v. We shall describe these 
excitations in terms of a field '?/'«(a;). This field describes particles and holes with momenta 
•p = fiv + /, where / -C /i. We will write I = Iq + l\\ + l± with l\\ = v{l ■ v) and /_l = / — /y. 
In order to take into account the entire Fermi surface we have to cover the Fermi surface 
with patches labeled by the local Fermi velocity, see Fig. |31 The number of such patches is 

~ (/i^/A^) where A_l <^ /i is the cutoff on the transverse momenta l±. 

Higher order terms are suppressed by powers of 1/ ^. As usual we have to consider all 
possible terms allowed by the symmetries of the underlying theory. At O ( l//i) we have 

-^^iDl^P.-ag^l^^^P^y (16) 

The coefficient of the first term is fixed by the dispersion relation of a fermion near the Fermi 
surface, Iq = l\\ + 1^/ (2/i) + . . .. The coefficient of the second term is most easily determined 
by writing the quark field in the microscopic theory as ip = tp^ + ip^ and then integrating 
out ip- at tree level. We find a = 1 + 0((7^), where the 0{g'^) terms arise from higher order 
perturbative corrections. At higher order in there is an infinite tower of operators of 
the form yU~""?/'^D^"^ (-u ■ Dy^'^^py with v = (1, —v) and n = 2ni + ^2 — 1. 
At 0(1/;U^) the effective theory contains four-fermion operators 



Vi r,r' 



The restriction + f 2 = f 3 + V4 allows two types of four-fermion operators. The first 
possibility is that both the incoming and outgoing fermion momenta are back-to-back. This 
corresponds to the BCS interaction 

^ = A E E vr'Rriv ■ ^) (^.r^-.) (^I'rvl,,) , (is) 

^ v,v' r,r' 

where v -v' = cos 6 is the scattering angle and Rf^'{v ■ v') is a set of orthogonal polynomials 
that we will specify below. The second possibility is that the final momenta are equal to 
the initial momenta up to a rotation around the axis defined by the sum of the incoming 
momenta. The relevant four-fermion operator is 

^ = A E E^r(0)^r(^^-^")(^.r^.')(^^rvl,), (19) 

where i), v' are the vectors obtained from v, v' by a rotation around Vtot = v + v' hy the angle 
0. In a system with short range interactions only the quantities -F)(0) are known as Fermi 
liquid parameters. The matrices F, F' describe the spin, color and flavor structure of the 
interaction. In the following we will focus on the spin structure. We can decompose 

(f) (f') = { (Hla^H^) {Hla^H^) , (Hla^aH^) (Hla^aH^) ] , (20) 

where H± = {1 ± v ■ (t)/2 are helicity projectors. In the spin zero sector there are two 

possible hehcity channels, (++) (++) and (++) — > ( ) together with their parity 

partners (-|- ^ — ). In the limit m ^ perturbative interactions only contribute to the 
helicity non-flip amplitude 

^ = AE^^^^K^- ^)(^.^2i^+^-.)(4^2i^+^L,) +(+---), (21) 

^^ v,v' 

where Pi{cos6) are Legendre polynomials. At 0{m? / ^'^) quark mass terms induce a non- 
zero helicity flip amplitude. The corresponding four-fermion operator was determined in 
[19^ . In Nf = 2 QCD instantons generate a helicity-flip amplitude which is suppressed by 
extra powers of {KqcdIij) l2Q]. In QCD with Nf = 3 flavors instantons produce a hel icity 
changing four-fermion interaction which is suppressed by both [Kqcd/iA ^^'^ ('^/a*) |22|- 
Even though helicity changing amplitudes are suppressed, they have important physical 
effects. For example, helicity flip amplitudes determine the masses of Goldstone bosons in 
the CFL and 2SC phase. 
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In the spin one sector there is only one hehcity channel (H — ) (H — ). The corresponding 
BCS interaction is 

C = ^^l^+-4?(t;- i/)(7/>,a2i/-ai/+^-.)(V'i;.a2i/-ai/+^l,,) + (+ ^ -), (22) 

where dii{cos6) is the reduced Wigner D-function. For details on the helicity amplitude 
formalism we refer the reader to 28|. Note that in the total helicity zero channel the reduced 
rotation matrix reduces to a Legendre polynomial, (IqqIcosO) = Pi{cos6), in agreement with 
equ. (PT|) . We will discuss the matching conditions for the the spin zero and spin one BCS 
helicity amplitudes VJ"*"^ and Vi^~ in Sect. IVIIII 

IV. HARD LOOPS 

There are two types of loop diagrams that arise in high density effective field theory, hard 
dense loops and soft dense loops. As an example of a hard dense loop we study the gluon 
two point function. At leading order in g and 1/yU we have 

n;t(P)^2.^iv,^Evv/|^ (,^_,^,(,;^^^_,^^^, . (23) 

where Ik = v- k. We note that taking the momentum of the external gluon to zero automat- 
ically selects forward scattering. We also observe that the gluon can interact with fermions 
of any Fermi velocity so that the polarization function involves a sum over all patches. After 
performing the fco integration we get 

where is the Fermi distribution function. We note that the integration is automatically 
restricted to small momenta. The integral over the transverse momenta l±, on the other 
hand, diverges quadratically with the cutoff A_|_. We observe, however, that the sum over 
patches and the integral over l± can be combined into an integral over the entire Fermi 
surface 

2tx^J (27r)2 271^ J An' ^ ^ 

This identification is consistent with our observation that the number of patches is n ~ 
/i^/(27r^A^). As a consequence of equ. the upper limit of the l± integral is effectively 
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fi, and not <^ /i- We will refer to loop integrals of this kind as hard loops. We find 

n;t(p)^2,.V'/^„„.„{i-^}. (26) 

This result agrees with equ. 0- We know however, that equ. (j^^ is not transverse and 
that the contribution from anti-particles, equ. Q is missing. In the effective theory, this 
contribution has to represented by a local counterterm. The required counterterm is [l5[ 

The appearance of this term is related to the fact that terms in the lagrangian of 
the form ^(w ■ D)"^^^, give non-vanishing tadpole contributions proportional to 

g2ni+n2 ^i-2ni~-n2 j^2ni _ ^-jn2 -yy-j-^jp}^ havc to bc absorbcd iuto couutertcrms. 

Putting everything together we find 

which agrees with the complete HDL result equ. (0). The gluonic three-point function can 
be computed in the same fashion. Using the rule equ. (^3]) we find 

TUP.^.r) ^ .,/-2™7f {^_|-^ - ^-l-^j . (29) 

which agrees with the HDL result equ. ()14j). We note that in the case of the three point 
function, as well as higher n-point functions, there is no leading order contribution from 
anti-particle states, and this is reflected by the absence of counterterms at leading order in 
the effective theory. 



V. SOFT LOOPS 

As an example of a soft loop contribution in the high density effective theory we study 
the fermion self energy. At leading order, we have 

m = I f \ ^ v,v.D,^ik), (30) 

J [ZttP Pq + ko - Ip+k 

where D^^{k) is the gluon propagator. Soft contributions to the quark self energy are 
dominated by nearly forward scattering. Note that this loop integral does not involve a sum 
over patches. Hard contributions to the fermion self energy are governed by the four-fermion 
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C = -m^ I^TrG,^-C^G^„ (31) 



operators equ. (fT^. We saw in the previous section that hard loops cause non-perturbative 
effects in gluon ra-point functions at the scale gfi. At energies below this scale we have to 
replace the free gluon term in the effective lagrangian by the generating functional for hard 
dense loops j29|, |30| 

^\ ^TrG^,^^_^^^G^, 

where the angular integral corresponds to an average over the direction of Pa = (IjP)- Note 
that the effective gluonic action is non-local. The effective fermion action in the high density 
effective theory for momenta below g^x remains local. The hard dense loop effective action 
equ. (|HTj) leads to the gluon propagator 

^-w-pij+A '''' 

where IIm and He are the transverse and longitudinal self energies in the HDL limit. The 
projection operators -P^^ are defined by 

Pjj = ^ij ~ kikj, Pqq = Pqi = 0, (33) 



P^u = -9,u + ^ - Pi- (34) 



In the regime |A;o| < \k\ < gfi the self energies can be approximated by H^; = 2m^ and 
Hm = i^m'^ko/\k\. We note that in this regime the transverse self energy is much smaller 
than the longitudinal one, 11^ < He- As a consequence the dominant part of the fermion 
self energy arises from transverse gluons. We have 

where Ip = v ■ p — fi and 1^ = v ■ k = Ix. To leading logarithmic accuracy we can ignore the 
difference between transverse and longitudinal cutoffs and set A_l = Ay = A. We can also 
set Ip ^ and (1 — x^) 1. We compute the integral by analytic continuation to euclidean 
space. Performing the integral over x we have 

S(p) = 2CEg' [ — [ arctan ( -i —. (36) 

The leading logarithmic term in the energy p4 can be extracted from 

. ,,^c. / |i / + p,) _ ^j^^^. (37) 

10 



Only the first term in t 



We get 



33 



le curly brackets gives a logarithmic contribution in the limit 0. 
35| 



" J + " ^''^ ' J ' ^ ^ 

which is accurate up to terms of 0{g^p4). Without calculating the 0{g^p4) term we cannot 
fix the scale inside the logarithm in equ. ()38|). 

Also note that the soft fermion self energy is of natural size. The soft gluon propagator 
scales as l/A^^^^, the soft fermion propagator scales as l/Agoft, and the loop integral gives 
AgQjj. As a result, the expected scaling is S ~ Ago ft which agrees with equ. (j38|l since p4 is 
a soft momentum. For comparison, the loop integral in the hard loop diagram discussed in 
section HVl scales as /i^A^^jj. 

Ward identities relate the soft fermion self energy to the soft quark-gluon vertex. The 
leading contribution to the abelian vertex function is given by 

r^f\puP2) = g'C^"'^ [ -^,v,v^v,D,,{k)- , (39) 

J (27r)4 (pi - A;)o - lp^_k {p2 - kjo - lp^_k 

with =Cf~ Ca/2 and Ca = N,. As in the case of the quark self energy the dominant 
contribution arises from the transverse part of the gluon propagator. Continuing to euclidean 
space we find 

dL r Pdl /-i l-x^ 



ia lPi,P2j 2vry (27r)2lr'^/| + /2 + |^2k 

— l)i — Ip^ — lxi{p2 — 1)4 — — Ix ^ ^ 
One can check that the term proportional to can be neglected. Performing the integral 
over X we get 

flog ( '^P' ' - ^^'^ + log ( '^''' - - ^^A ] m 

Brown et al. observed that there are two distinct kinematic regimes for the vertex function, 
depending on whether {pi — ^2)4 is bigger or smaller than {Ip^ — Ip^) ^J. In the limit 

Q'Pi - ^2) < (Pi - ^2)4 we get 

(42) 
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where (^1)4 — (^2)4 — P4- This integral is clearly identical to equ. (jHTj) for the derivative of 
the fermion self energy. We get 

hm hm T^f\p^,P2) = ^l^log - • (43) 

In the opposite limit (/p^ — Ip^) ^ (pi —^2)4 the delta function is not present and the vertex 
function does not have a logarithmic divergence. 

We saw that in the hard dense loop approximation the abelian and non-abelian vertex 
corrections are identical up to a color factor. The sum of the color factors of the two diagrams 
is C('^f') + = {Cf - Ca/2) + Ca/2 = Cp which shows that the quark-gluon vertex in 

the HDL approximation has the same color structure as the quark self energy. The same 
result is also obtained for the leading IR divergence in the soft dense loop limit. The only 
difference is that following the arguments given above we have to take into account both 
the HDL gluon self energy and the HDL gluon three-point function when evaluating the 
non-abelian vertex function, see Fig. infc. Indeed, the free gluon three-point function does 
not contribute to the IR divergence of the quark-gluon vertex in the SDL limit. We have 

r^'HPi,P2) = g'C-^^'^J 0-^v,v,D,,ik)T^,,,,{p2-puk,p,-p2-k) 

■D,>,{pi -p2 - k)— 1 . (44) 

Again the interesting limit is (/p^ ~ lp2) ^ ipi ~ ^2)4- The calculation is greatly simplified 
by the observation that in this limit equ. gives 

^4fiu{P2 - Pi,k,pi - P2 - k) ~ i — — {Il^^{k + {p2 - pi)) - Ilf,^{k)} 

{P2-Pi)4 

This relation is also a direct consequence of the Ward identity for the HDL gluon three-point 
function. Using equ. (j^5|) and continuing to euclidean space we get the following expression 
for the non-abelian vertex function in the limit (/p^ — ipi) ^ {pi ~ ^2)4 

This integral can be evaluated using the same strategy as for the abelian vertex. First we 
observe that to leading logarithmic accuracy we can replace (1 — x^) \. The integral over 
X then gives 

1 2i ( I \ 

Ti ^ ;~ = T arctan . (47) 

tliA + PAj-lx I \U+PaJ 
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In the limit p4 — the integral over I4 is dominated by the discontinuity of the inverse tan 
function and 

hm hm ri-'')(Pi,P2) = ^^^log - . (48) 

The non-abelian vertex function in the limit {Ip^ ~ ^p2) ^ {Pi ~ ^2)4 can be computed 
using the same methods. We find that there is no logarithmic enhancement near the Fermi 
surface. We conclude that the abelian and non-abelian soft dense loop vertex functions have 
the same logarithmic singularity and that the sum of the two contributions is proportional 

to C^'^'') + = Cp. 



VI. COLOR SUPERCONDUCTIVITY 



Soft dense loop corrections to the fermion self energy and the quark-gluon vertex function 
become comparable to the free propagator and the free vertex at the scale ~ /x exp(— l/gf^). 
This implies that at this scale soft dense loops have to be resummed. Physically, this 
resummation corresponds to the study of non- Fermi liquid effects in dense quark matter . 
However, before non-Fermi liquid effects become important the quark-quark interaction in 
the BCS channel becomes singular. The scale of superfluidity is ~ /iexp(— !/(/) j^^. This 
scale is smaller than the HDL scale but larger than the SDL scale. This implies that we 
have to resum quark-quark scattering with HDL dressed gluon propagators but that SDL 
corrections to the quark self energy and the quark-gluon vertex are small and can be treated 
perturbatively. 

The resummation of the quark-quark scattering amplitude in the BCS channel leads to 
the formation of a non- zero gap in the single particle spectrum. We can take this effect into 
account in the high density effective theory by including a tree level gap term 

C = /\R^i{v ■ k)ij_^a2T^^ + h.c.. (49) 

Here, F is any of the helicity structures introduced in Sect. IIIIl i?f [x) is the corresponding 
angular factor and A is a unit vector. The magnitude of the gap is determined variationally, 
by requiring the free energy to be stationary order by order in perturbation theory. We 
shall see that the gap varies with the energy or the residual momentum on a scale set by 
the gap itself. As a result the energy dependence is relevant and we cannot replace the gap 
by its value on the Fermi surface. In order to do perturbation theory in the presence of 
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a gap term we will use the Nambu-Gorkov method and introduce a two component field 
\E' = {ipy, The inverse propagator for the field \1/ is 

S-U,)J ^W^'^f"--^)^), (50) 

\A*{p)Rf{v- A)T^ v-p J 

where = {1,—v). The variational principle for the gap A gives the Dyson-Schwinger 
equation 

2g^ f d\ A(g4) ^ ( \ 

^^^'^ = — J j2^^ ql + P^ + Aigy ^''''^^^P - 

where the factor in front of the integral is the color factor corresponding to the color anti- 
symmetric [3] channel. It is sufficient to solve this equation to leading logarithmic accuracy. 
We shall see in the next section that corrections to the gap can be computed perturbatively, 
without solving the gap equation to higher accuracy. To leading logarithmic accuracy the 
gap equation is dominated by the IR divergence in the magnetic gluon propagator. This IR 
divergence is independent of the helicity and angular momentum channel. We have 

18^' ^0 Jql + A{q^y \\pl-qlV'v 



The leading logarithmic behavior is independent of the ratio of the cutoffs and we can set 
Aj| = = A. We introduce the dimensionless variables variables x = log(2A/(g4 + e^)) and 
y = log(2A/ (p4 + ep) where Cg = {ql + A^q^))^^"^. In terms of dimensionless variables the gap 
equation is given by 

^^y^ = T^ dxA{x)K{x,y), (53) 

loTT^ JO 

where Xq = log(2A/Ao) and K{x,y) is the kernel of the integral equation. At leading order 
we can use the approximation K{x,y) = min{x,y) |3S]. We can perform an additional 
rescaling x = XqX, y = x^y. Since the leading order kernel is homogeneous in x, y we can 
write the gap equation as an eigenvalue equation 

A{y)=xl-^ dxA{x)K{x,y), (54) 

where the gap function is subject to the boundary conditions A(0) = and A'(l) = 0. Son 
observed that equ. (fH^ is equivalent to the differential equation (2^ 
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which has the solutions 

An{x) = Anfl sin (^ 3^27]- ^°'"^) ' ^ ^"^^ ~^ V^' ^^^^ 

The physical solution corresponds to ri = which gives the largest gap, Aq = 
2A exp(— 37r^/(v^5'))- Solutions with 7^ have smaller gaps and are not global minima of 
the free energy. 



VII. PERTURBATIVE CORRECTIONS TO THE GAP PARAMETER 

The complete set of solutions of the leading order gap equation can be used to set up 
a perturbative scheme for computing corrections to the gap function. This perturbative 
scheme is similar to the one used by Brown et al. in order to compute corrections to the 
critical temperature j^. We first observe that the functions A„(x) form an orthogonal set 
of solutions of the gap equation 

/ dx An{x)Am{x) = 5m,n- (57) 

Jo 

We now consider the gap equation with the kernel K{x, y)+6K{x, y) where K{x^ y) contains 
the leading IR divergence and 5K{x, y) is a perturbative correction. We have 

^(y) = ^07^ A{x)K{x,y) + XQ-—T dx A{x)6K{xox,xoy). (58) 

We write the gap function A{x) and the eigenvalue perturbative expansion in 6K 

A{x) = A(o)(x) + A(^)(x) + A(2)(x) + . . . , (59) 
xo = x^o^ + xi'^ + xj,'^ + . . . . (60) 

Using the orthogonality of the unperturbed solutions A„(x) we can derive expressions for A^*) 
and Xq\ These expressions are very similar to ordinary Rayleigh-Schroedinger perturbation 
theory. At first order in perturbation theory we get 



^(1) _ 1 (A0)\^ 9 



and 



IStt^ Jo Jo 



7^-^4°^t|-^ ['dx ['dyAS'\x)6KixoX,Xoy)Af\y), 

1 _ ( 1 V ISvr^ JO JO 

\.2k+lJ 



2 rl rl 

dx / dy A''Q\x)6K{xoX,Xoy)A''Q\y), (61) 



(62) 
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where c^^^ are the components of /\^^\x) in the basis of the unperturbed gap functions, 
/S.'^'^\x) = J2Ck^^^f\x). Equ. shows that the first order correction to the pairing gap 
can be computed using the unperturbed gap function. 

The simplest example for a correction to the kernel is is a contact term 

5K{xox,xoy) = log(6). (63) 

In Sect. lVliTl we will show that the coefficient b is determined by four-fermion operators in the 
effective theory. The unperturbed gap is given by A = Aexp(-4°^) = Aexp(-37rV(y2^)). 
Using equ. (jHT|) and equ. we get 

A = Ae-(4°'+4'^+-) = 6Ae"^. (64) 

We conclude that contact terms modify the prefactor of the gap. We also notice that because 
of the orthogonality of the gap functions a contact term does not modify the shape of the 
gap function. 

We can also study the effect of the quark self energy Equ. ()38|1 implies a wave 

function renormalization 

ZW-l-f^>o.(^). (65) 

The corresponding correction to the kernel of the gap equation is 

9^ - 

5K{xoX, xov) = -T-^ {xqx) K{xoX, Xoy), (66) 

where we have used Cp = 4/3 for Nc = 3 and K(x,y) = mm{x,y) is the leading order 
kernel. Using equ. we find 

where we have used the fact that Ao^-* (x) is an eigenfunction of the unperturbed kernel. The 
result equ. (jU7j) corresponds to a reduction of the gap by a factor exp(— (4 + 7r^)/8). We can 
also see that the fermion self energy correction leads to an admixture of higher harmonics 
of the gap function. However, these admixtures are small, c^.^"* = 0{g), for all values of k. 

There is a slight subtlety with regard to the result equ. (jHZI)- In section |V| we computed 
the wave function renormalization in the normal phase. We observed, however, that the 
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scale of non-Fermi liquid effects, E ~ fi exp{—l / g"^) , is exponentially small as compared 
to the scale where pairing sets in, E ~ fiexp{—l/g). This implies that also the normal 
component of the fermion self energy should be computed with the gap taken into account. 
It is easy to see that in this case equ. (j38p is modified to 



Pairing removes the infrared divergence in the wave function renormalization. However, in 
the superfiuid phase the fermion self energy is only modified for very small energies < A 
whereas the correction to the gap is dominated by the regime > A. As a consequence 
the result equ. ()67|) is not changed. Finally, we have to consider the role of soft dense 
loop vertex corrections given in equ. PHj) and ()48|) . The infrared logarithm in the vertex 
correction only appears in the regime \lp — lq\ ^ \p4 — q4,\. Since the gap equation determines 
the anomalous self energy on the quasi-particle mass shell, p^ ~ (/^ + A^)^/^, this implies 
that Ig <^ g4. However, this condition eliminates the BCS logarithm in the gap equation 
fl51|) . As a consequence, the infrared logarithm in the quark-gluon vertex does not modify 
the eigenvalue Xq at 0(1) in the coupling constant. Vertex corrections to the gap equation 
were first considered in [39] but the arguments given in that work are not correct. A more 
detailed study can be found in 3^. 

VIII. MATCHING CONTACT TERMS 

Consider the effect of the 0(l//i^) contact term equ. (PT|) on the pairing gap. The 
contribution of this term to the kernel of the gap equation involves a sum over patches, see 
Fig.lHl As a consequence the correction to the kernel, g'^6K{x,y) = 6V(^^ , is not suppressed 
by 1/yU^. In the previous section we showed that 0{g^) constants in the kernel contribute to 
the eigenvalue at 0(1). This implies that we have to determine to coefficients of 0{g'^/fi'^) 
contact terms in the high density effective theory in order to determine the eigenvalue of 
the gap equation to 0(1). 

This can be achieved by matching the quark-quark scattering amplitude in the BCS 
channel. The tree level scattering amplitude in the spin zero color anti-triplet channel is 
given by 




(68) 




(69) 
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At leading order in the effective theory this amphtude is represented by 

1 1 



HDET 



(0) 



3 



+ 



2/^2(1 _ cos ^) +ni5 2/^2(1- COS ^) + nM 



l-cose<A2 /(2ai2) 



+ 



-P;(cOSI 



(70) 



where the colhnear term is cut off at = A^. The matching condition requires that the 
partial wave amplitudes corresponding to equ. and (ffn|) are the same. Since the collinear 
term in the high density effective theory is regulated by a UV cutoff the counterterm depends 
on the cutoff, too. The matching condition is simplest for A^^ = 2/i2. The s-wave term is 
given by 



3 4 



dx 



-i(l-x) 



+ 



h(l-x) 



0, 



(71) 



(1 - a;) + ns/(2/i2) (1 - x) + nM/(2/i2) 

up to corrections of 0{g^). The cutoff dependence of Vf^'^ is controlled by the renormaliza- 
tion group equation 



A^ 



d 



dAl 



3_ 

3 



(72) 



This equation imp 



We find 



39, 



4G 



ies that the kernel in the gap equation is independent of the cutoff Aj^. 



K{pi, Qa) = log 



A, 



+ log 



/ 5127r' 



7^ 



5A 



Nf 



(73) 



Using the leading order result equ. (fHUI). the perturbative corrections equ. (j64|67p . and the 
value of the counterterm given in equ. (f7H|l we get the standard result for the gap in the 
2SC phase 

A = 5127r^(2/Ar^)^/V^7"'e-^e"^. (74) 

Note that at this order we did not have to fix the dependence of the counterterms on the 
longitudinal cutoff Ay. 

We saw that Vo^^(2/i2) vanishes to leading order in perturbation theory. This is not the 
case for higher partial wave terms and operators with non-zero spin. For example, the total 
angular momentum one counter term in the total helicity zero channel is 



'3 

>2 2 



I)] -_2 ^ 



X 



(75) 
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The tree level matrix element in the helicity one channel is proportional to (1 + cos^)/2 for 
both electric and magnetic gluon exchanges. The corresponding counterterm is 

4?(^)[(5 + f) + (Uf)]-2 



dx ■ 



1 — X 



9V 
2 3 



t;^^ (76) 



where we have used d^ii{x) = (1 +x)/2. These results imply that the total angular momen- 
tum one gaps are given by Af^ = exp(—6)An~^ and = exp(— 9/2)A([^ where A^^ is 
the s-wave gap given in equ. (fTijl (3lL l42l |43|. 



IX. SUMMARY 



In this work we studied several problems related to the use of effective field theories in 
QCD at large baryon density. We showed that the power counting in is complicated by 
the appearance of two scales inside loop integrals. Hard dense loops involve the large scale 

and lead to phenomena such as screening and damping at the scale gfi. Soft loops only 
involve small scales and lead to superfluidity and non- Fermi liquid behavior at exponentially 
small scales. We also showed that contact terms in the effective lagrangian are suppressed by 
powers of but they get enhanced by hard loops. As a result contact terms contribute 
to the pre-exponent of the pairing gap at 0(1). We performed the necessary matching 
calculation to determine four-fermion operators in the effective theory. 

There are many problems that remain to be studied. We would like to understand not 
only how to count powers of E/fi but also logarithms. This is essential for understanding the 
structure of the expansion of the gap and other quantities in the coupling constant. Another 
issue is the choice of regularization scheme. In this work we have used a momentum space 
cutoff throughout. There are, however, some HDET calculations that have been performed 
using dimensional regularization. It is clearly important to understand the power counting 
in both schemes. 

We would also like to have a better understanding of gauge invariance. The high den- 
sity effective theory has the advantage that the pre-exponent of the gap is determined by 
counterterms that are matched against on-shell scattering amplitudes and that are therefore 
manifestly gauge invariant. This result agrees with explicit calculations in a generalized 
Coulomb gauge ^J] but not with calculations in a general covariant gauge [4^. A formal 
argument that the full quasi-particle propagator is gauge invariant was recently presented in 
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4^- It is well known that the hard dense loop effective action equ. ()3ip is gauge invariant. 
It is not clear whether there is a similar statement for soft dense loops. Brown et al. showed 
that the soft quark self energy and quark gluon vertex function satisfy BRST identities 0| . 
We have also checked that the IR divergences in the soft wave function renormalization and 
vertex are independent of the gauge parameter in a general Coulomb gauge. On the other 
hand, Hsu et al. have argued that one can find a non-local gauge in which there is no wave 
function renormalization 

Hong and Hsu have argued that the high density effective theory can be formulated 
non-perturbatively, using a lattice regulator, and that the leading order theory has a pos- 
itive euclidean measure j47|. Such a formulation might be useful in order to study certain 
non-perturbative questions and to provide guidance in the search for algorithms capable of 
simulating gauge theories at non-zero baryon density. Our results suggest, however, that 
there are some difficulties with this proposal. The leading order effective theory is defined on 
a single patch, or, if a gap term is added, on two patches corresponding to Fermi velocities 
±v. In this theory the ground state is likely a coherent particle-hole state of the type sug- 
gested by Deryagin, Grigoriev and Rubakov The reason is that the forward scattering 
amplitude in the particle-hole channel is larger than the one in the particle-particle channel. 
This state is disfavored in the full theory, because the corresponding four-fermion operator 
is not enhanced by hard dense loops. This means that operators that are superficially sub- 
leading are important in selecting the correct ground state. Physically, this is related to the 
fact that the particle-hole state can only become coherent over the entire Fermi surface if 
the surface is nested, but not if the Fermi surface is spherical. Another problem is how to 
correctly incorporate the boundary conditions for the effective theory. The lagrangian of the 
leading order theory is equal to the non-relativistic QCD (NRQCD) lagrangian. However, 
the physics of the high density theory and NRQCD are very different. For example, there 
is no screening in the gluon two-point function in NRQCD. This issue is also related to the 
necessity to add a two-gluon counterterm to the effective theory. In perturbation theory this 
difference is encoded in different ie prescriptions for the fermion propagator. It would be 
interesting to understand how the boundary conditions are realized on a euclidean lattice. 
Some of these issues are studied in 
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c) 




FIG. 1: Hard dense loop contributions to the the gluon and quark self energy. Fig a) shows the 
particle-hole contribution to the gluon self energy, b) shows the anti-particle contribution, and c) 
shows the quark self energy. 

a) b) 




FIG. 2: Hard dense loop contributions to the gluon three and four function. 
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FIG. 3: Hard dense loop contribution to the quark-gluon vertex function. Figure a) shows the 
abeHan vertex and figure b) the non-abehan vertex. 




FIG. 4: High density effective field theory description of excitations near the Fermi surface. 
The effective theory is defined on patches labeled by the local Fermi velocity v. Momenta are 
decomposed with respect to v, p = ^xv + + 



24 



high energy scale 

E gH 

screening and damping 

E ~ /iexp(-l/^) 
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Fermi surface 



FIG. 5: Scales of the high density effective field theory, 
a) 




FIG. 6: Soft dense loop contribution to the quark self energy (Fig. a)) and the quark-gluon vertex 
function (Figs, b) and c)). The solid squares denote the HDL gluon self energy and three-point 
function. 
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FIG. 7: Leading order gap equation for the superconducting gap in the high density effective 
theory. The soUd square denotes the HDL self energy and the open square is the anomalous quark 
self energy. 



a) b) 




> A > □ < i < 

FIG. 8: Corrections to the gap equation. Fig. a) shows the soft dense loop quark self energy 
correction, Fig. b) local four-fermion interactions, and Fig. c) shows the soft dense loop vertex 
correction. The solid triangles denote the SDL self energy and vertex correction and the solid 
circle is the four-fermion interaction in the effective theory. 
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